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On  the  Influence  of  Internal  Density  Variations  on  the 
Linear  Stability  Characteristics  of  Planar  Shear  Layers 

A.  Hajesfandiari  and  D.J.  Forliti1 
Department  of  Mechanical  and  Aerospace  Engineering 
The  State  University  of  New  York  at  Buffalo 
Buffalo,  NY 

Planar  shear  layers  are  known  to  be  sensitive  to  density  variation  within  the  vorticity  region.  The 
current  study  will  expand  on  the  available  literature  on  the  stability  characteristics  of  shear  layers 
with  density  variation  through  consideration  of  new  regimes  in  tenns  of  the  velocity  and  density 
profiles.  Both  spatial  and  spatio-temporal  linear  analyses  were  considered.  In  addition  to  the 
maximum  density  change  within  the  shear  layer,  the  effect  of  density  to  vorticity  thickness  ratio, 
and  transverse  shifting  of  the  density  profile  relative  to  the  velocity  profile  were  considered.  For 
a  shear  layer  experiencing  combustion,  these  effects  relate  to  the  relative  position  and  thickness 
of  the  flame  in  the  shear  layer.  Thick  density  profiles  have  a  uniformly  stabilizing  affect  that 
imposes  essentially  a  low-pass  attenuator  on  the  spatial  stability.  The  shear  layer  is  destabilized 
through  the  presence  of  an  intense  yet  thin  density  reduction  that  overlaps  with  the  peak  vorticity 
profile,  a  condition  that  may  be  present  near  the  shear  layer  origination  point.  The  velocity  ratio 
at  which  the  flow  transitions  between  absolute  and  convective  instability  is  dependent  on  the 
magnitude  of  the  density  reduction  in  the  shear  layer  as  well  as  the  relative  density  profile 
thickness  and  transverse  location.  The  results  generally  indicate  that  the  instability  of  the  planar 
shear  layer  is  strongly  influenced  by  the  characteristics  of  the  density-weighted  vorticity  within 
the  shear  layer. 


I.  INTRODUCTION 

Linear  stability  theory  has  been  found  to  be  a  powerful  tool  for  predicting  parametric  effects  of 
shear  flows  including  the  effects  of  velocity  ratio1'3,  density  ratio4'6,  heat  release7'9,  and 
compressibility10'12.  The  results  of  such  linear  studies  provide  a  theoretical  prediction  of  initial 
instability  frequency  or  wavelength,  the  instability  growth  rate,  and  the  convective/absolute 
nature  of  the  instability,  knowledge  that  is  critical  for  understanding  and  controlling  shear  flows. 
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A  number  of  linear  stability  analyses  have  been  applied  to  reacting  shear  flows.  Shin  and 
Ferziger  and  Mahalingam,  Cantwell,  and  Ferziger  conducted  detailed  linear  stability  studies  of 
reacting  shear  flows  including  perturbations  in  heat  release.  These  studies  found  that  the  effects 
of  heat  release  were  accurately  replicated  when  only  accounting  for  the  base  velocity  and  density 
profiles.  In  further  support  of  this  approximation.  Day,  Reynolds,  and  Mansour11  suggest  that 
the  mean  density-weighted  vorticity  profile  is  a  strong  indicator  of  instability  for  compressible 
reacting  shear  flows.  It  is  apparent  that  this  approximation  is  likely  to  breakdown  for  reacting 
flows  near  quenching  boundaries  . 

The  primary  objective  of  the  current  research  was  to  investigate  the  affects  of  heat  release  in  the 
shear  layer  region  on  the  convective/absolute  instability  transition.  Experimental  studies  indicate 
that  crossing  from  convective  into  absolute  instability  is  associated  with  the  emergence  of  a  self- 
excited  state  that  experiences  rapid  spatial  growth,  entrainment,  and  energetically-enhanced 
turbulence14'16.  Although  heat  release  in  general  has  a  stabilizing  effect  on  the  shear  layer,  and 
reduces  entrainment17'22,  there  is  an  opportunity  to  use  linear  theory  to  identify  reacting  shear 
layer  configurations  that  are  absolutely  unstable.  A  reacting  shear  layer  in  an  absolute  instability 
regime  will  likely  have  enhanced  entrainment  characteristics,  providing  an  opportunity  to 
increase  combustion  intensity.  Additionally,  absolutely  unstable  flows,  if  properly  tuned,  have 
potential  use  in  combustion  control  applications  to  avoid  thenno-acoustic  instabilities.  This 
analysis  will  also  promote  understanding  of  the  role  of  flame  position  and  flame  thickness 
relative  to  the  shear  layer  on  the  instability  characteristics.  In  particular,  the  flame  thickness 
relative  to  the  shear  layer  may  be  very  different  in  the  flame  stabilization  region  compared  to  the 
fully-developed  turbulent  regime  located  further  downstream. 

We  denote  the  two-dimensional  base  flow  velocity  components,  density  and  pressure  by  U  ,  V  , 
p  and  P  ,  and  the  corresponding  quantities  for  the  disturbance  by  u  ,v ->P  and  p  .  The 
problem  can  be  simplified  by  the  assumption  of  parallel  base  flow,  i.e.  U  -  U(y )  , 
V  =  0,  p  =  p{y)  and  p  =  const .  The  results  of  Day  et  al.11  indicate  that  at  low  convective  Mach 
numbers,  the  shear  layer  with  heat  release  will  be  dominated  by  two-dimensional  instabilities, 
thus  three-dimensional  instability  was  not  considered. 

The  solution  to  the  linear  disturbance  equations  can  be  decomposed  into  normal  modes  .  The 
transverse  perturbation  velocity  term,  for  example,  becomes  of  the  fonn 
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v'{x,y,t)  =  v(y)  exp  [i(ax-(Ot)\ , 


(1) 


where  CL  is  the  wavenumber  and  CO  is  the  frequency.  The  complex  wave  speed  is  defined  as 

co 

c  =  cr+ict  =  —  . 

a 

The  substitution  of  the  normal  modes  into  the  governing  equations  and  reducing  the  system  of 
equations  to  a  single  equation  in  tenns  v  yields 
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For  the  planar  shear  layer,  asy  — >  ±oo ,  equation  (2)  simplifies  to  the  form 


'Vv 

dy2 


a 2  v 


=  0. 


(3) 


which  goes  to  zero  exponentially  as  y  — >  ±  co . 

Equation  (2)  together  with  the  boundary  conditions  (Eq.  3)  and  base  profiles  fonn  the  stability 
model  for  two-dimensional  incompressible  parallel  shear  layer  with  density  variation.  The 
system  of  equations  is  solved  using  the  shooting  method. 


II.  Shear  Layers  with  Internal  Density  Variation 

,_Ul-U2 

The  shear  layer  velocity  ratio  is  defined  as  lI  —  — — — ,  where  U]  and  U2  are  the  high  and 

Ul  +  U  2 

low-speed  freestream  velocities,  respectively.  The  velocity  profile  is  in  the  form  of  a  hyperbolic 
tangent  shape,  the  most  common  model  for  analyzing  free  shear  layers24 

U(y)=(MII2)tmh(2y/Sa)+U ,  (4) 

where  A  U  is  the  velocity  difference  across  the  layer,  U  is  the  average  velocity,  and  8a  is  the 
vorticity  thickness  defined  as 

A  U 


5(0  U'  ’ 


(5) 
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where  U'm ax  is  the  maximum  slope  of  the  mean  velocity. 

In  order  to  study  the  stability  of  shear  layers  with  heat  release,  computations  are  developed  for 
shear  layers  with  a  density  reduction  within  the  shear  layer  region.  In  addition  to  considering  the 
intensity  of  the  density  change,  the  relative  thickness  and  transverse  location  of  the  density 
reduction  were  also  considered.  The  density  profde  used  for  this  purpose  is 

p(y)  =  1  +  A  tanh (By  -0.1 )  — A  tanh(f?  y  +  O.l)  (6) 

A  -  1-5 

where  A  -  — -  and  s  is  the  ratio  of  the  minimum/maximum  density  to  the  free  stream 

2  tanh  (0.1)  y 

density, 


s  = 


Po_ 

Poo 


(7) 


The  thickness  of  the  density  profile  is  adjusted  by  varying  the  parameter  B.  A  measure  used  to 
assess  the  thickness  of  the  density  profile  is  the  density  thickness,  8  .  This 

quantity  is  defined  as 


Sp=2 


_Ap_ 

P  max 


(8) 


where  A p  is  the  magnitude  of  the  density  change.  The  variable  r  is  used  to  represent  the 
relative  density  thickness  and  is  defined  as 


r 


(9) 


The  assumed  profile  shapes  are  used  to  facilitate  a  general  parametric  study  of  the  influence  of 
heat  release  on  a  planar  shear  layer.  Real  planar  flames  will  exhibit  deviations  from  the  assumed 
shape,  including  the  presence  of  asymmetry  in  the  velocity  and/or  density  profiles.  In  general, 
the  results  of  linear  stability  with  assumed  profile  shapes  have  been  predictive  of  the  results  of 
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real  experimental  flows,  as  demonstrated  by  Strykowski  and  coworkers14'16’  25.  It  is  thus 
anticipated  that  the  results  will  capture  the  general  trends  of  the  affects  of  heat  release. 

III.  Influence  of  the  Density  Ratio  s 

The  spatial  analysis  has  been  performed  for  the  base  velocity  and  density  profiles  shown  in  Figs. 
1  and  2,  respectively.  In  Fig.  2,  density  profiles  are  plotted  versus  the  vorticity  magnitude  profile 


Figure  1.  Base  velocity  profile  U (y)  =  1.5  tanh  (y)  +1.5 


Figure  2.  Base  density  profile  versus  the  vorticity  magnitude  profile  for  internal  density  variation  for  r  =3.0 
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for  a  relative  density  thickness  r  =  3.0  .  The  dimensionless  spatial  growth  rate,  -  cet  ,  versus 


the  dimensionless  frequency  co  ,  is  illustrated  in  Fig.  3  for  different  density  ratios.  The 
dimensionless  frequency  and  growth  rate  are  defined  as 


co  - 


A  U 


(10) 


and 


=  «  A 


(ll) 


An  increase  in  density,  which  is  destabilizing,  is  also  shown  for  reference.  The  figure  indicates 


Figure  3.  Spatial  growth  rate  for  internal  density  variation  for  r  =3.0 
the  stabilizing  effect  of  reducing  the  density  ratio.  The  low-frequency  side  of  the  growth  rate 
profile  is  nominally  unchanged  by  the  presence  of  heat  release  until  the  peak  is  approached.  It  is 
also  apparent  that  the  high-frequency  side  of  the  growth  rate  profile  experiences  a  constant 
reduction  in  amplification  rate  as  indicated  by  the  parallel  curves.  The  stabilizing  effect  of  heat 
release  is  essentially  a  low-pass  attenuator  with  a  cutoff  frequency  that  depends  on  the  density 
ratio,  and  an  asymptotically  constant  reduction  magnitude  at  high  frequency;  the  damping  effect 
does  not  grow  with  increase  in  frequency  at  high  frequency.  For  density  ratios  above  unity,  the 
density  variation  acts  as  a  low-pass  amplifier  of  instability. 
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The  same  analysis  was  done  for  a  case  with  the  same  velocity  profile  but  a  density  profile 
corresponding  to  a  thickness  ratio  r  =  1.0.  This  base  flow  may  be  expected  to  be  associated  with 
the  region  near  the  splitter-plate  trailing  edge  of  a  reacting  shear  layer  under  conditions  where 
the  combined  boundary  layer  thicknesses  of  the  two  streams  is  larger  than  the  diffusion  flame 
thickness.  Although  r  is  unity,  Fig.  4  shows  that  the  density  variation  is  indeed  thinner  than  the 

10  f  -  S  =  0.25 

-  S  =  0.5 

-  S  =  1.0 

- S  =  2.0 

-  _  .  Vorticity 


_l _ i _ i _ i _ i _ _ i _ i _ i _ i _ I _ i _ i _ i _ i _ l _ i 

0.5  1  1.5  2 

p  ,  dU/dy 

Figure  4.  Base  density  profile  versus  the  vorticity  magnitude  profile  for  internal  density  variation  for  r  =  1 .0 

vorticity  variation.  Figure  5  clearly  demonstrates  a  similar  stabilizing  trend  with  decreasing 
density  ratio  except  for  s  =  0.25.  By  decreasing  the  density  ratio  from  0.5  to  0.25,  the  peak 
spatial  growth  rate  increases  dramatically.  The  case  for  s  =  0.5  shows  some  increase  in  instability 
at  low  frequencies  that  was  not  seen  in  Fig.  3.  It  is  clear  from  Fig.  5  that  the  low-pass  attenuator 
affect  changes  for  the  thinner  density  profile  for  sufficiently  low  density  ratio.  This  transition 
correlates  to  changes  in  the  shape  of  the  density-weighted  vorticity  profile,  which  will  be 
described  in  the  next  section. 

IV.  Influence  of  the  Thickness  Ratio 

The  nonmonotonic  behavior  illustrated  in  Fig.  5  clearly  relates  to  the  relative  thickness  of  the 
density  profile,  since  it  is  not  observed  in  the  results  for  the  thicker  density  profile  (Fig.  3). 
Therefore,  the  next  step  was  to  investigate  the  affect  of  the  relative  thickness  on  the  stability  of 
the  flow.  Figure  6  shows  the  amplification  rate  of  the  most  unstable  mode  at  s  =  0.5  and  0.25  for 
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Figure  5.  Spatial  growth  rate  for  internal  density  variation  for  V  =  1.0 

different  values  of  the  thickness  ratio  r.  It  is  very  clear  from  the  figure  that  the  influence  of  heat 
release  on  the  stability  of  the  shear  layer  is  complex.  At  sufficiently  high  values  of  r,  the  heat 
release  is  clearly  stabilizing  (growth  rates  are  smaller  than  the  uniform  density  case).  There  does 
exist  a  range  of  r  at  low  values  where  the  heat  release  is  actually  destabilizing.  Further  insight 
into  this  is  gained  through  consideration  of  the  density-weighted  vorticity  profile. 

As  suggested  by  Day  et  ah,  the  density-weighted  vorticity  plays  a  key  role  in  the  instability 
characteristics  of  shear  layers.11  At  a  fixed  (but  less  than  unity)  density  ratio,  the  density- 
weighted  vorticity  profile  may  have  two  different  shapes.  At  large  r,  the  density  weighted- 
vorticity  profile  exhibits  a  single  centered  peak.  In  this  regime,  the  affects  of  heat  release  are 
uniformly  stabilizing.  At  small  r,  the  density-weighted  vorticity  profile  will  exhibit  two  peaks 
located  symmetrically  away  from  the  center  of  the  shear  layer.  These  two  regimes  are  indicated 
on  Fig.  6.  In  the  two-peak  regime,  the  heat  release  can  have  a  destabilizing  effect.  This  regime 
is  also  associated  with  low-frequency  amplification  as  shown  in  Fig.  5;  the  low-pass  attenuator 
behavior  is  associated  with  the  single  peak  regime. 

Figure  7  shows  the  peak  magnitude  of  the  density-weighted  vorticity  as  a  function  of  r  and  s.  In 
the  single  peak  regime,  the  peak  amplitude  of  the  density-weighted  vorticity  does  not  depend  on 
the  thickness  ratio.  In  the  two-peak  regime,  a  reduction  in  thickness  ratio  actually  results  in 
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Figure  6.  Spatial  amplification  rate  of  the  most  unstable  mode  as  a  function  of  r 


Figure  7.  Peak  density-weighted  vorticity  as  a  function  ofr 

higher  amplitudes  of  the  peak  density-weighted  vorticity,  which  according  to  Day  et  al.  suggests 
enhanced  instability.  This  is  consistent  with  the  trend  shown  in  Figs.  5  and  6.  The  results  of 
Figs.  6  and  7  further  suggest  the  utility  of  the  density-weighted  vorticity  profile  as  an  indicator  of 
instability.  Although  the  peak  density-weighted  vorticity  for  s  =  0.25  is  always  lower  than  that 
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Figure  8.  Density  profile  shift  with  respect  to  the  velocity  profile  for  r  =  2.0 

for  s  =  0.5,  very  thin  density  layers  for  s  =  0.25  are  much  more  unstable  than  for  s  =  0.5,  as 
shown  in  Fig.  6.  For  these  very  thin  layers,  the  smaller  density  ratio  s  will  tend  to  push  the 
somewhat  lower  amplitude  peak  further  away  from  the  center  of  the  shear  layer.  Thus  the 
overall  instability  is  dependent  both  on  the  peak  magnitude  in  the  density-weighted  vorticity  as 
well  as  the  transverse  location  of  the  peak(s).  Further  transverse  displacement  of  the  vorticity 
peak  appears  to  result  in  enhanced  instability. 

V.  Influence  of  Shifting  the  Density  Profile 

In  addition  to  the  influence  of  density  layer  thickness,  the  affect  of  transverse  shifting  of  the 
density  profile  was  explored.  The  analysis  consists  of  shifting  the  density  profile  while  keeping 
the  velocity  profile  fixed.  For  a  shear  layer  experiencing  combustion,  the  shift  relates  to  the 
position  of  the  flame  in  the  shear  layer  and  is  illustrated  by  Fig.  8  for  r  =  2.0.  Physically,  the 
shifting  of  the  heat  release  location  in  the  shear  layer  is  related  to  the  lateral  location  of  the 
stoichiometric  mixture  fraction,  which  is  governed  by  the  reactant  concentrations  in  the 
freestreams  and  the  entrainment  ratio.  Shifts  towards  both  the  high-  and  low-speed  streams  were 
considered.  The  magnitude  of  the  displacement  is  expressed  in  relation  to  the  vorticity  thickness, 
8m .  Figure  9  shows  the  amplification  rate  for  different  offsets  of  the  density  profile  for  r  =  3.0 
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Figure  9.  Spatial  growth  rate  for  internal  density  variation  with  shifted  density  profile,  r  =  3.0 ,  S  =  0.25 

and  s  =  0.25.  This  case  is  in  the  single  density-weighted  vorticity  peak  regime.  As  the  density 
profile  is  moved  towards  the  upper  or  lower  stream,  the  flow  becomes  more  unstable.  It  should 
be  noted  here  that  shifting  the  density  profile  in  this  single  peak  regime  of  the  density-weighted 
vorticity  results  in  a  higher  peak  value  (slightly  skewed  away  from  the  shifting  direction),  thus 
suggesting  enhanced  instability.  It  is  apparent  that  shifting  the  density-weighted  vorticity 

towards  the  low-speed  side  has  a  more  dramatic  destabilizing  influence.  Results  for  shifting 
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other  thickness  ratio  s  cases  may  be  found  in  Hajesfandiari  . 

VI.  Spatio-temporal  Instability 

The  idea  of  absolute  and  convective  instability  was  first  introduced  by  Briggs  and  Bers“  in  the 
study  of  plasma  instabilities.  The  absolute-convective  nature  of  the  instability  is  identified  when 
there  exists  at  least  one  point,  where  da>/da=  0  ,  in  the  complex  CC  plane.  This  point  is  a 
singularity  point  in  the  wavenumber  space  and  a  saddle  point  in  the  frequency  space. 

In  one  of  the  initial  applications  of  this  concept  towards  fluid  mechanics,  Huerre  and 
Monkewitz  applied  the  spatio-temporal  analysis  to  hyperbolic  tangent  shear  layers  with 
uniform  density.  They  solved  the  stability  equation  for  different  values  of  velocity  ratio  to 
explore  the  affect  of  this  parameter  on  the  absolute-convective  transition  of  the  instabilities. 
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Based  on  their  findings,  the  zero  group  velocity  mode  with  zero  5  i.e.  oj-  -  0 ,  occurs  when 
A  =  1.315  .  This  critical  velocity  ratio  corresponds  to  U2/Ul  =-0.136  .  This  implies  that 
absolute  instability  may  start  to  emerge  in  the  shear  layer  in  the  presence  of  a  sufficiently  high 
counterflow.  Strykowski  and  Niccum15  experimentally  studied  the  counterflowing  mixing  layer 
for  different  velocity  ratios.  The  critical  value  of  the  velocity  ratio  obtained  by  them  was 
A  =  1.32  which  agrees  strongly  with  linear  theory. 

In  the  present  work,  exploring  the  role  of  the  internal  density  variation  on  the  absolute- 
convective  instability  transition  of  free  shear  layers  is  of  interest.  For  this  purpose,  the  spatio- 
temporal  analysis  is  applied  to  different  flow  profiles  in  order  to  find  out  at  what  value  of 
velocity  ratio  the  transition  from  convective  to  absolute  instability  occurs  in  the  presence  of 
internal  density  variations.  Additionally,  the  effect  of  the  density  thickness  and  transverse  offset 
on  the  transition  velocity  ratio  was  also  studied. 

Figure  10  shows  the  critical  velocity  ratio  at  which  the  instability  characteristics  switch  from 


Figure  10.  Critical  velocity  ratio  as  a  function  of  density  ratio 


absolute  to  convective.  The  results  are  represented  for  two  density  profiles  that  are  different  in 
thickness.  For  thick  density  layers,  absolute  instability  is  inhibited  (i.e.  requires  higher  levels  of 
counterflow)  except  at  very  low  density  ratios,  while  the  trend  for  the  thinner  density  layer 
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appears  to  promote  absolute  instability  with  decreasing  density  ratio.  The  trend  of  a  thinner  layer 
promoting  absolute  instability  is  in  agreement  with  the  findings  of  Reynal  et  al.  for  low-density 
planar  jets.  For  density  ratios  below  0.25,  the  critical  velocity  ratio  for  r  =  1  is  smaller  than 
unity,  suggesting  absolute  instability  without  the  requirement  of  counterflow.  This  scenario  has 
potential  applicability  to  the  flame  anchoring  region,  either  attached  to  the  trailing  edge  or  lifted, 
where  the  flame  and/or  flammable  region  may  be  thinner  than  the  vorticity  layer  of  the  velocity 
profile.  The  vorticity  and  density  thickness  in  the  downstream  region  are  more  likely  to  be  the 
same  due  to  turbulent  transport  (visually  this  would  be  associated  with  r  ~  2),  thus  the  flame  may 
have  limited  impact  on  the  transition  velocity  ratio.  In  general,  the  results  suggest  that  local 
regions  of  absolute  instability  may  be  promoted  or  abated  depending  on  the  relative  flame 
thickness,  but  the  global  flame  which  may  include  large  regions  of  thickness  ratio  near  two  will 
have  approximately  the  same  tendency  to  be  absolutely  unstable  (thus  globally  unstable)  as  a 
nonreacting  shear  layer. 

The  effect  of  transverse  shifting  of  the  density  profile  on  the  absolute-convective  transition  of  the 
instability  was  also  explored.  The  value  of  the  critical  velocity  ratio  with  shifting  of  the  density 
profile  to  the  high  and  low-speed  streams  is  illustrated  in  Fig.  11  for  r  =  3.0.  When  the  density 


A/8 


Figure  11.  Critical  velocity  ratio  as  a  function  of  transverse  offset  when  r  =3.0 


profile  is  shifted  to  the  low-speed  side,  absolute  instability  is  less  prevalent  in  that  higher  levels 
of  counterflow  are  required.  This  stabilizing  effect  of  shifting  the  density  profile  to  the  low- 
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speed  side  is  consistent  with  the  results  of  Reynal  et  al.  and  Srinivasan,  Hallberg,  and 

t  i 

Strykowski  for  low  density  jets.  Divergent  from  the  trend  shown  by  Reynal  et  ah,  shifting  the 

density  profile  to  the  high-speed  side  results  in  enhancement  of  absolute  instability.  Similar 

26 

effects  are  observed  for  different  thickness  ratios  and  can  be  found  in  Hajesfandiari . 

X.  SUMMARY 

In  general,  the  results  indicate  that  the  influence  of  a  layer  of  heat  release  on  a  planar  shear  layer 
can  have  a  wide  range  of  affects  on  the  hydrodynamic  instability.  When  the  density  layer  is 
relatively  thick  compared  to  the  vorticity  thickness,  a  reduction  in  density  in  the  shear  layer  has  a 
stabilizing  effect  in  terms  of  spatial  growth.  If  the  layer  is  shifted  towards  one  of  the 
freestreams,  it  appears  that  the  stabilizing  effect  is  somewhat  degraded,  and  the  flow  can  be 
highly  unstable  when  the  heat  release  is  shifted  to  the  high-speed  side.  The  affect  of  heat  release 
is  essentially  a  low-pass  attenuator  that  degrades  amplification  above  a  cutoff  frequency  that 
depends  on  density  ratio. 

For  the  case  of  a  relatively  thin  density  profile,  the  reduction  of  the  density  to  a  sufficient  level 
results  in  an  eventual  destabilization  of  the  flow,  with  growth  rates  that  can  exceed  the  uniform 
density  case.  This  destabilization  effect  degrades  if  the  layer  is  shifted  away  from  the  center  of 
the  shear  layer. 

The  above-stated  conclusions  further  support  the  utility  of  the  density-weighted  vorticity  as  an 
indicator  of  instability.  The  peak  value  within  the  density-weighted  vorticity  profile  is  sensitive 
to  the  thickness  ratio  when  two  peaks  are  present  in  the  profile,  and  the  instability  appears  to  also 
be  very  sensitive  to  thickness  ratio  r  in  this  regime.  In  the  range  of  profile  characteristics  where 
a  single  peak  is  present,  the  peak  density-weighted  vorticity  is  constant  with  changing  thickness 
ratio,  and  the  instability  is  very  weakly  dependent  on  the  thickness  ratio.  In  this  regime,  the 
thickness  ratio  causes  very  slight  changes  in  the  profile  shape  of  the  density-weighted  vorticity. 
One  should  be  cautioned  to  not  place  too  much  emphasis  on  the  density-weighted  vorticity.  The 
instability  characteristics  will  be  very  much  dependent  on  the  freestream  conditions  as  well  as 
the  behavior  of  the  density-weighted  vorticity  profile.  One  case  in  point,  the  instability  of  a 
hyperbolic  tangent  velocity  profile  with  a  freestream  density  difference  will  be  highly  dependent 
on  whether  the  low-density  side  is  on  the  high  or  low  velocity  side  of  the  shear  layer.  These 
cases  would  have  nominally  the  same  density-weighted  vorticity  profile.  Nonetheless,  processes 
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internal  to  the  shear  layer  that  alter  the  distribution  of  the  density-weighted  vorticity  can  have 
dramatic  consequences  on  the  stability  without  changes  in  the  freestream  conditions. 

The  results  provide  insight  on  the  instability  of  different  spatial  regions  of  reacting  shear  layers. 
In  general,  the  turbulent  region  of  the  shear  layer  will  be  associated  with  a  density  profile  that  is 
of  similar  thickness  to  the  velocity  profile,  due  to  the  governing  nature  of  the  turbulent  transport. 
Exceptions  would  include  the  case  where  the  stoichiometric  mixture  fraction  is  dramatically 
different  from  the  entrainment  ratio  of  the  shear  layer,  which  would  lead  to  a  flame  that  resides 
on  the  outside  (or  inside)  edge  of  the  shear  layer,  or  exhibiting  asymmetry.17'19  If  the  flame  is 
located  nominally  within  the  shear  layer  and  has  similar  thickness,  the  flame  tends  to  stabilize 
the  shear  layer,  although  transition  to  absolute  instability  occurs  at  a  similar  level  of  counterflow. 
This  provides  encouragement  that  global  instability  may  be  achievable  in  reacting  shear  layers 
with  realizable  counterflow  levels.  It  is  apparent  that  different  behaviors,  including  higher 
growth  rates  and  lower  transitional  velocity  ratios,  are  present  when  the  density  profile  is  thinner 
than  the  velocity-based  shear  layer  profile.  This  affect  for  thin  flames  is  very  sensitive  to  the 
transverse  location  of  the  flame.  The  results  in  general  support  the  need  for  further  exploration 
of  employing  global  instabilities  in  reacting  flows  to  enhance  mixing  and  combustion  processes. 
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